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Propagation  of  Radius-Tailored  Laser  Pulses  Over 
Extended  Distances  in  a  Uniform  Plasma 

I.  Introduction 

The  physics  of  intense  laser  pulse  propagation  in  an  underdense  plasma  is  a  problem 
of  recent  interest  for  a  variety  of  applications/*^  including  radiation  generation  (laser 
harmonics,  x-ray  sources)  and  advanced  accelerators  (laser  wake-field  accelerator,  plasma 
beat-wave  accelerator).  Some  applications,  particularly  the  laser  wake-field  accelerator 
(LWFA),  require  laser  propagation  over  several  laser  diffiraction  lengths  at  ultra-high  laser 
intensities.  In  the  LWFA,*“*  a  short  (tl  <  1  ps),  high  power  (P  >  1  TW)  laser  pulse 
propagates  in  plasma  to  generate  a  large  amplitude  {Eg  >  1  GV/m)  wake  field,  which  can 
trap  and  accelerate  a  trailing  electron  bunch. 

Ultra-high-intensity  laser  pulses  require  a  tight  focus  and  are  characterized  by  an 
intensity  high  enough  that  relativistic  effects  become  significant  (oo  ^  1)  and  a  vacuum 
diffraction  length,  or  Rayleigh  length  Zjt,  that  is  is  rdatively  short  {Zr  <  1  cm).  Here, 
Oo  —  eAo/m«e*  is  the  peak  amplitude  of  the  normalized  vector  potential  of  the  laser 
pulse,  Zr  =  irTg/A,  ro  is  the  laser  spot  size  at  focus,  A  is  the  laser  wavelength,  m,  is 
the  electron  mass,  and  a  Gausrian  radial  profile  has  been  assumed.  In  terms  of  the  peak 
laser  intenrity  It  oo  ^  8.5  x  10~^"A[/tm]p/^[W/cm’]  for  linear  polarization,  which  will  be 
assumed  throughout  this  paper. 

Computer  simulations  of  such  laser  pulses  can  be  problematical  because  of  the  dis¬ 
parate  time  and  length  scales  involved.  TirpicaHy,  X  <.  L  C  Za,  where  the  laser  pulse 
length  L  =  cri  b  defined  as  the  fuU-width-at-half-maximum  of  the  laser  intensity.  A 
numerical  code,  LEM  (Laser-ElectroMagnedc),  has  been  developed  to  simulate  laser  prop¬ 
agation  in  an  underdense  plasma  for  arbitrarily  high  values  of  oq.  It  b  based  on  the 
cold-fluid  model  equations  f^ven  in  Re&.  7  and  8.  In  developing  thb  code,  the  model  equa¬ 
tions  £ar  the  plasma  response  to  a  fpven  laser  field  vrere  recast  in  a  numerically  tractable 
fixrm.  In  thb  paper,  the  numerical  equations  are  presented  and,  as  an  example,  the  code 
b  applied  to  the  propagation  of  radius^taUored  laser  pulses.  .Theoretical  model  equations 
govemiag  the  choice  of  radius  tailoring  are  abo  fpven. 
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Radius-tailored  laser  pulses,  which  were  discussed  briefly  in  Ref.  7,  have  a  unique  prop¬ 
erty  in  that  they  can  propagate  over  extended  distances  (>  IOZa)  in  a  uniform  underdense 
plasma.  While  previous  simulations  of  such  pulses^  were  performed  in  the  context  of  the 
paraxial  approximation,  which  pves  a  constant  laser  group  velocity  Vg  =  c,  the  present 
study  shows  that,  for  cases  of  interest,  variations  in  the  group  velocity  of  the  laser  pulse 
over  the  length  of  the  pulse  can  have  a  significant  effect.  In  the  present  study,  theoreti¬ 
cal  model  equations  governing  the  chmce  of  the  anal  pulse  profile  for  radius  tailoring  are 
discussed  and  parameter  requirements  and  tolerainces  are  ^ven  for  such  pulses  and  for  the 
injection  of  these  pulses  into  a  plasma.  Also,  the  construction  of  a  radius-tailored  laser 
pulse  from  from  a  series  ultra-short  Gaussian  pulses  is  considered. 
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n.  Model  Equations 

These  simulations  were  based  on  the  laser-plasma  fluid  model  described  in  Re&.  7 
and  8,  which  utilises  (r,^  =  z  —  et^r  =  <)  coordinates.  The  laser  pulse  moves  in  the 
positive  z  direction  such  that  the  &ont  of  the  laser  pulse  remains  near  ^  =  0.  The  physical 
region  of  interest  extends  from  C  =  0,  where  the  plasma  is  unperturbed,  to  ^  <  0.  The 
model,  which  takes  advantage  of  the  separation  between  the  fast  (A)  and  slow  (ZR,Ap) 
time  and  space  scales,  is  valid  when  Zr  L,  Zr  Xp,  A  <  ro,  and  A  <  A^,  where 
Xp  =  Tirtfuip  is  the  plasma  wavelength,  =  (49re^no/me)^/^  is  the  electron  plasma 
ficequency,  no  =  n^®^(r  =  0)  and  n^®^(r)  is  the  initial  electron  dennty  profile  (in  the 
simulations  shown  below,  an  initially  uniform  plasma  radial  density  profile  is  assumed). 
Laser  pulse  evolution  is  described  by  the  wave  equation 


A 

e  dr 


(1) 


1V>  indude  variations  in  the  laser  group  vdodty,  the  term  is  retained  in  Eq. 

(1),  in  contrast  to  Ref.  7.  In  Eq.  (1),  =  eA//m.e*  is  normalised  vector  potential  of 

the  laser  pulse,  •/  is  the  slowly  varying  amplitude  (a/  s  a/eiq>(»fcoC)/2  +  c.e,  where  c.c. 
denotes  complex  coqjngate),  kp  —  w^/c,  p,  as  n«/7«no>  n,  is  the  plasma  density,  7,  is  the 
relativistic  factor  of  the  plasma,  and  the  sulMcript  s  denotes  a  slowly-varying  component. 
The  plasma  response  to  a  fpven  laser  fidd  ay  is  pven  by*^’® 


and 


Vl«.  =  t^.u.  - 

(2o) 

'i  +  ^  {f-P-  -  />'*’) . 

(2i) 

^(a.-a,)  =  V(T,-#,), 

(2c) 

2(1 -h^.)  • 

(2d) 

tdiete  the  CSonlomb  gauge  has  been  used  (V  *  a«  »  0)  and  ip,  =  a,,,.  In  Eqs.  (2a- 

d),  ac  SB  eA«/me^  and  ss  are  the  normalised  vector  and  scalar  potentiab. 
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respectively,  u«  =  Pt/mc  is  the  normalized  momentum,  =  p,{(  =  0),  and  the  plasma 
ions  are  assumed  to  be  immobile. 

In  the  axisymmetric  case,  Eqs.  (2a-d),  along  with  V  •  a«  =  0,  can  be  combined  to  yield 
a  single  equation^**  of  the  form  This  equation  can  be  written  as 

^  =  (*J/>.-Vl)u.„  +  A(VxUx..),  (3.) 

where  u,.,  =  7,  -  -  1, 

and 

A  form  of  equation  (3a)  is  solved  numerically  to  obtain  the  plasma  response.  Both  the 
derivation  of  Eqs.  (3a-c)  and  the  details  of  the  numerical  solution  are  discussed  in  the 
Appendix. 

The  LEM  code  has  been  used  to  simuiate  the  effects  of  several  instabilities  of  recent 
interest,  such  as  the  forward  Raman  scatterix^  (FRS)  instability*  and  the  self-modulation 
instalnfity.***'*  Certain  instabilities,  however,  cannot  be  adequately  described  wing  the 
formulation  fpven  here.  In  particular,  the  growth  large-angle  Raman  scattering,  also  called 
Raman  nde-scatter  (RSS),^*’^^  is  suppressed  by  the  numerical  scheme,  which  includes  a 
radial  smoothing  algorithm.  This  smoothing  algorithm  is  required  to  suppress  a  numerical 
instalnlity  and  is  discussed  in  the  Appendix. 

It  is  reasonable  to  neglect  the  effect  of  RSS  on  the  simulation  results  given  below 
because  of  the  short  pulse  duration  {h^L  s'  12)  and  because  the  intensity  varies  over  the 
length  of  the  pulse.  SpedficaUy,  Ref.  10  indicates  that  the  growth  of  RSS  over  the  lexigth 
of  the  pulse  is  minimal  for  short  pulses  with  i^L  <  10  —  20.  Furthermore,  additional 
simulatioxis  of  radius-tailored  pulses,  with  kpL  >  10,  showed  that  the  growth  of  the  FRS 
and  sdf-xxiodulation  instabilities  were  stroxi^  suppressed  by  the  radius  tailoring.  This 
suggests  that  radius  taifoxing  may  also  suppress  large-angle  Raman  scattering. 
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in.  Radius-Tailored  Laser  Pukes 

A  radius-tailored  laser  pulse  is  a  pulse  configuration  which  is  designed  to  to  propagate 
over  a  long  distance  by  the  mechanism  of  relativistic  optical  gtiiding.^^*^’  Relativistic  opti¬ 
cal  guiding}  which  generally  occurs  when  the  laser  power  Pq  —  21 .5(00^0 /A)^  GW  exceeds 
the  critical  power  Pc  2;  17(Ap/A)*  GW,  where  Po/Pe  =  (ibproao)*/32,  is  problematical  be¬ 
cause  the  leading  portion  of  a  laser  pulse  (|^|  <  Ap)  tends  to  dif&ract  as  if  in  vacuum^*^  and 
because  long  pulses  {L  >  Ap)  tend  to  undergo  self-modulation.*'^”  With  radius  tailoring, 
the  User  spot  sise  at  focus  is  tapered  firom  a  large  value  at  the  pulse  front  to  a  small  value 
over  the  length  of  the  laser  pulse  while  keeping  P/Pe  1  throughout  the  pulse.  This 
condition  is  met  when  a£(^)rx,(0  —  0.9Ap,  where  is  the  laser  spot  size  and  ol  =  |d/|. 
Radius  tailoring  is  derigned  to  reUtivistically  guide  the  high  intensity  portion  of  the  laser 
pulse  while  allowing  the  leading  portion  of  the  pulse  to  simply  diffract.  The  tapered  radius 
produces  a  locally  Urge  diffraction  length  at  the  pulse  front,  Za,o(C  =  0)  >  Zr,o(C  = 
where  Zj|,o  =  irr£(^,r  =  0)/A.  Hence,  the  pulse  can  remain  focused  over  Urge  distances 
cr  >  Zji,o(C  ~  Furthermore,  self*moduUtion  is  suppressed  for  L  <  2Ap  because  the 
increaang  User  intensity  over  the  length  of  the  pulse  does  not  excite  a  strong  wake  field 
within  the  pulse.  Note  that  while  the  spot  vise  at  focus  varies  along  the  length  of  the  pulse, 
each  point  along  the  puke  must  reach  its  focus  simultaneously. 

The  sdection  of  a  puke  profile  k  motivated  by  considering  the  envelope  equation  of 
ReC  9,  which  describes  the  evolution  of  the  User  spot  sise  r^((,r)  of  a  puke  with  an 
iqqproadmately  Gaussiaa  radial  profile  in  the  limits  a*  <  1,  >  1.  Assuming  an 

initially  uniform  plasma,  the  envelope  equation  can  be  written  as 

^  ^  {1  -  =  ».  (4) 

where  M  ss  rj;((,r)/ro,  rg  k  now  defined  to  be  a  constant  equal  to  the  TnitiiTniim  spot 
rise  at  focus  (corresponding  to  the  peak  value  of  |a/|  —  og)  and  Zr  =  srJ/A.  The 
second  tom  in  Eq.  (4)  represents  vacuum  difi^action  whereas  the  third  term  represents 
Use  focuring/diffiractive  effects  of  the  plasma. 

Initially,  for  cr  <  |Zji^//|,  R  m  R(t  a=  0)[1  +  where  the  inverse  of 
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(5) 


the  eSective  diffraction  length  is 


PiOlPc 


HO  +  RHC 


')]’ }  ■ 


hi  this  notation,  k  —  l/Zj{,o  corresponds  to  vacuum  diffraction,  k  =  0  indicates  a 
‘‘matched”  ^-slice  (focusing  and  diffraction  forces  are  balanced)  and  <  0  corresponds 
to  net  focusing.  The  chmce  of  initial  pulse  profile,  specified  as  P(^)  and  ai,((,T  =  0),  is 
motivated  by  the  desire  to  have  k  vary  smoothly  from  k  =  l/Z/{,o(C  =  0)  at  the  leading 
edge  of  the.  puke  to  k  ~  0  at  the  tail  of  the  pulse. 

To  illustrate,  consider  four  radius-tailored  puke  configurations,  each  at  constant 
power:  (a)  ai  =  oo  ^(—(ir/4Xp),  0  <  — C  <  2Xp,  P  =  P^;  (b)  Case  (a)  with  P  =  1.2Pe; 
(c)  oj  =  — Cao/2Ay,  0  <  — ^  <  2Ay,  P  =  P*;  (d)  Case  (c)  with  P  =  1.14Pc.  Plots  of 
(kZjr)^  versus  ^  are  shown  for  each  case  in  Fig.  1.  Here,  (kZr)^  =  ^R/^A,e//(C)  ^ 

measure  of  the  initial  rates  of  diffraction/focusing  throughout  the  puke  relative  to  the 
vacuum  diffraction  rate  associated  with  the  back  of  the  puke.  For  comparison  {kZr)^  for 
a  rectangular  puke  with  ai,  =  oo,  0  <  — C  <  2Ap,  and  P  =  Pc,  k  shown  in  Fig.  2.  Figures 
1  and  2  iUiutrate  key  phenomena  that  have  been  observed  in  LEM  simulations  of  tailored 
and  non-taikxred  laser  pukes. 

Hie  most  significant  difiaction  in  Fig.  1  k  observed  at  the  peak  of  (kZr)^  near  the 
beam  head  (|C|  <  Ay/2),  where  the  User  k  everywhere  slowly  diffracting.  For  example, 
the  peak  diffraction  in  case  (a)  occurs  at  C  —  — 0.65Ay,  where  {kZr)*  =  0.015.  Thk 
corresponds  to  Zr,«//  =  8.2Zjt.  The  peak  diffraction  near  the  beam  head  can  be  reduced 
by  increadng  the  laser  power  in  thk  repon.  In  case  (b),  the  peak  diffraction  corresponds 
to  Zji^ff  =  10.5Zjt.  Thk  suggests  that  a  carefully  deagned  tailored  puke  might  have 
peak  power  somewhat  greater  than  Pe  at  and  near  the  puke  head,  tapering  to  P  Pe  at 
the  tafi  of  the  puke. 

Fm  constant  P/Pey  there  k  a  trade-off  between  di&action  near  the  head  and  ovor- 
foenring  in  the  body  of  the  pulse.  Simulations  show  that,  at  ro  =  the  allowable  range 
k  1  <  P/Pe  <  1*5,  with  i^timum  results  at  P/Pe  1.2.  While  a  detailed  study  has  not 
been  perCmned,  it  k  dear  that  thk  tolerance  will  vary  with  ro.  In  particular,  for  ro  >  Ay 
and  P  >  Pe,  have  verified  that  the  rdativktically-focused  tail  portion  of  the 
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pulse  is  far  from  radial  equilibrium.  For  P  >  Pc,  a  laser  pulse  in  a  transverse  equilibrium 
state^’’^*  has  a  characteristic  radius  <  Xp. 

The  tendency  of  non-tailored  laser  pulses  to  erode^'^  is  illustrated  in  Fig.  2,  which 
indicates  di&action  over  a  region  of  length  Xp/4  near  the  head  of  the  rectangular  pulse 
(^R,c//  =  Zr  &t  the  pulse  head).  Figure  2  also  shows  regions  of  focusing  and  defocusing 
which  lead  to  the  sdf-modulation  instability  of  Refs.  6-10.  These  periodic  regions  are  due 
to  the  plasma  wave  generated  by  the  rise  time  £rt«e/c  of  the  pulse  (for  a  square  pulse, 
^  1).  Figure  1  shows  that  radius— tailored  pulses  also  have  this  tendency,  but  to  a 
lesser  degree  (the  characteristic  defocusing  lengths  are  smaller  by  approximately  an  order 
of  magnitude).  Simulations  show  that  for  L  ^  Xp  self-modulation  can  be  a  significant 
effect,  even  for  radius-tailored  laser  pulses. 
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rv.  Simulations 


As  an  example  of  radius-tailored  laser  pulse  propagation,  consider  a  laser  pulse  with 
radius  tailoring  over  the  region  0  <  — C  <  2Ap,  where  P  =  1.14Pc  and  a\  =  -Cao/2Ap  as 
in  case  (c)  above,  and  decreasing  power  over  the  region  2Ap  <  —(  <  5Ap/2,  where  r^,  =  rp 
and  o|_  =  — 2oq(^  +  2Ap)/Ap.  With  oo  =  0.96  and  plasma  density  no  =  5.5  x  10^^  cm"’, 
we  have  ro  =  Ap  =  45  /xm,  Pg  =  35TW,  Po  =  40TW,  and  Zr  =  457rAp  =  0.64cm. 

The  simulation  begins  at  r  =  0  with  the  laser  pulse  outside  the  plasma,  one  Zr  away 
firom  the  focal  point.  The  transition  from  vacuum  to  full  plasma  density  occurs  over  a 
distance  of  0.3Za  such  that  full  density  is  reached  at  the  focal  point.  Simulations  show 
that  with  a  longer  transition  refpon,  the  laser  over-focuses  and  the  effectiveness  of  radius 
tailoring  is  reduced.  The  simulation  continues  until  er  =  16Za. 

Figure  3  shows  plots  of  the  normalised  laser  intensity  |a/p  at  er  =  Zr,  as  the  pulse 
enters  the  plasma  with  a  minimum  spot  size  r^.min  =  Ap,  and  at  er  =  16Za,  where  the 
minimum  spot  sise  is  =  0.59Ap.  In  this  code,  the  ri  is  the  radius  enclosing  86.5% 

of  the  laser  power  (for  a  Gaussian  pulse,  a  ~  ). 

Several  interesting  effects  have  been  observed  in  this  and  in  other  LEM  simulations. 
Firstly,  the  laser  pulse  deviates  from  its  initial  Gaussian  profile,  becoming  highly  peaked 
on  aads.  As  a  result,  the  quantity  ax,rj;  does  not  remain  constant.  Secondly,  because 
Vg  <  c,  there  is  a  noticeable  degree  of  slippage  in  the  ^-position  of  the  pulse  versus  time. 
Thirdly,  because  of  nonlinear  effects,  Vg  varies  with  the  laser  intensity.’’^*  The  tendency 
of  the  laser  to  ‘‘steepen”  and  develop  a  peaked  profile  versus  C  can  be  attributed  to  this 
effect.  In  frict,  the  range  of  validity  of  these  simulations  was  limited  by  this  tendency: 
when  ai\dai,ld^\~^  ^  A,  Eqs.  (2a)-(2d)  no  longer  hold.  Fourthly,  a  large  amplitude  wake 
fidd,  with  Em  ~  30  GV/m,  is  generated  behind  the  laser  pulse. 

Note  that  the  results  stated  above  and  in  the  previous  section  indicate  requirements 
on  plasma  unifbrmity  and  on  the  sharpness  of  the  vacuum-plasma  boundary.  Because 
Pc  ~  1/ns,  the  requirement  that  1  <  P/Pg  ^  1*5  indicates  that  a  variation  in  no  as  large 
as  ±15%  is  acceptable  for  a  radius-tailored  pulse  with  uniform  power.  This  was  verified 
to  some  dqpree  in  simulations  of  radius-tailored  pulses  in  which  the  plasma  density  varied 
by  ±10%  with  2-4  oscillations  per  Zjt*  These  denmty  variations  produced  no  discernible 
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change  in  the  results.  When  the  laser  power  is  nonuniform,  as  in  the  cases  below,  this 
requirement  becomes  somewhat  more  strict.  The  requirement,  stated  above,  that  the 
vacuum-plasma  interface  have  a  width  smaller  than  O.ZZr  should  not  pose  a  significant 
difficulty,  since  Zr  Xp  for  typical  tailored-pulse  parameters. 

While  the  ideal  laser  pulse  configuration  of  Fig.  3  may  be  difficult  to  produce  in  the 
laboratory,  it  may  be  possible  to  construct  an  approximate  radius-tailored  pulse  through 
the  superposition  of  several  ultra-short  pulses,  each  with  a  different  spot  size.  A  somewhat 
idealized  configuration  of  5  overlapping  laser  pulses  is  shown  in  Fig.  4.  Each  pulse  has  a 
Gaussian  radial  profile  and  an  axial  profile  ^ven  by  ai  sin(— ^x/2X),  0  <  — C  <  2X,  with 
X  =  17  ftm  (L/c  =  57  fs).  This  configuration  is  idealized  in  the  sense  that  the  overlapping 
pulses  are  phase-locked,  resulting  in  constructive  interference,  and  the  separation  between 
pulses  is  23  fjan.  such  that  the  variation  in  power  along  the  pulse  is  minimized  with  peak 
power  Pp^ak  Po  =  40  TW,  as  in  Fig.  3,  and  (-averaged  power  (P)  ~  Pc.  This  simulation 
is  otherwise  identical  to  that  shown  in  Fig.  3.  As  in  that  case,  a  laige  amplitude  wake  field 
Eg  ^  25  GV/m  is  generated  behind  the  lastf  pulse. 

A  less  ideal  configuration  of  overlapping  short  pulses  is  illustrated  in  Fig.  5.  Here,  the 
idealization  of  phase-locking  is  removed,  with  each  pulse  having  a  randomly-chosen  phase 
shift.  Since  tailored  pulse  simulations  at  ro  =  Xp  show  successful  guiding  for  1  <  P/Pc  < 
1.5,  the  pulse  spacing  has  been  decreased  to  20  /im.  Thus,  in  the  case  of  perfect  constructive 
interference,  Ppgmh  —  l<5Po«  where  again  Po  =  40  TW.  Because  this  configuration  contains 
a  randmn  phase-shift  for  each  pulse,  the  results  from  an  ensemble  of  simulations  must  be 
conndered.  For  practical  reasons,  a  statistically  snudl  ensemble  of  20  simulations  was  used. 
Of  these,  five  were  deemed  ‘‘successful”,  i.e.,  a  wake  with  Eg  >  20  GV/m  persisted  over 
a  propagation  range  greater  than  10  Zr.  Results  from  one  such  case  are  ^ven  in  Figs. 
5  and  6,  with  fig.  5  showing  the  laser  intennty  plotted  at  cr  =  Za  and  at  cr  =  IIZa, 
and  Fig.  6  showing  the  laser  power  plotted  versus  (  at  r  =  0.  For  comparison,  P  versus 
(  is  also  shown  for  the  ideal  phase-locked  case  of  Fig.  3.  The  results  of  the  remaining  15 
simnlatiims  were  characterized  by  the  axial  break-up  of  the  laser  pulse.  In  many  of  these 
cases,  one  or  more  of  the  resulting  “beamlets”  propagated  greater  than  10  Za,  but  with  a 
reduced  wake  ampfitude  (Eg  <  15  GV/m). 
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Clearly,  experimental  confirmation  of  the  results  of  Figs.  4  and  5  would  be  difficult. 
While  it  is  beyond  the  scope  of  this  paper  to  provide  an  experimental  design,  it  seems 
reasonable  to  speculate  that  in  such  an  experiment  a  long  laser  pulse  would  be  split  into 
five  separate  pulses  with  separate  optical  paths.  The  optics  would  have  to  be  timed  in 
such  a  way  that  the  pulses  arrive  at  focus  simultaneously.  Only  small  errors  in  the  relative 
timing  and  focal  positions  (<!  20  /xm)  would  be  allowed.  It  is  not  clear  whether  or  not 
such  an  experiment  could  be  performed  using  present  technology. 


V.  Conclusions 


The  numerical  fluid  code  LEM,  described  above  and  in  the  Appendix,  allows  simula¬ 
tion  of  intense  laser  pulses  over  extended  distances  (cr  >>  Zn).  Because  the  code  computes 
the  evolution  of  the  laser  envelope  on  the  plasma  time-scale,  these  simulations  are  signifi¬ 
cantly  faster  than  those  carried  out  by  more  conventional  means,  such  as  the  particle-in-ceU 
method.  Simulations  of  radius-tailored  laser  pulses  propagating  in  a  uniform  underdense 
plasma  show  that  such  pulses  can  remain  focused  over  distances  greater  than  IOZr  in 
accordance  with  £q.  (5),  which  was  derived  firom  the  laser-plasma  envelope  equation  of 
Ref.  9.  These  simulations  show  that  large  amplitude  wake  fields  (Eg  >  1  GV/m)  can  be 
generated  behind  the  radius-tailored  pulse  provided  that  the  intensity  terminates  in  less 
than  a  plasma  wavelength.  Simulations  also  show  that  a  radius-tailored  laser  pulse  can  be 
constructed  firom  a  series  of  overlapping  ultra-short  (Z/c  50/s)  Gaussian  ptilses. 
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Appendix:  Numerical  Solution  of  the  Model  Equations 
The  model  equations  are  given  as  Eqs.  (2a-d)  in  the  main  text.  These  equations  make 
use  of  the  quasi-static  approximation,  in  which  it  is  assumed  that  the  scale  length  over 
which  the  laser  field  evolves  is  long  compared  to  the  laser  pulse  length  {Zr  ^  L).  In  this 
limit  derivatives  with  respect  to  r  are  dropped  relative  to  derivatives  with  respect  to  ^  in 
the  plasma  fluid  equations,  but  not  in  the  wave  equation.  In  addition,  Eqs.  (1)  and  (2a-d) 
have  been  averaged  over  the  fast  (A)  laser  oscillations.  The  resulting  simulation  code, 
LEM  (Laser-EHectroMagnetic),  describes  the  laser  envelope  d/  and  the  plasma  response 
on  a  spaitial  grid  that  resolves  Ap,L  and  ro,  with  a  time  step  that  resolves  the  evolution  of 
the  laser  intensity  {Zr/c)  smd  the  phase  osdllations,  due  to  plasma  dispersion,  of  the  real 
and  itn aginary  parts  of  d/,  which  have  a  characteristic  period  kX^firc. 

Equations  (2a-d)  are  numerically  intractable  in  the  sense  that  it  is  not  clear  which 
equation  should  be  solved  first,  second,  etc.,  or  if  an  iterative  method  of  solving  the  equa¬ 
tions  can  be  expected  to  converge  to  a  meaningful  solution.  This  numerical  quagmire  can 
be  avoided  by  recasting  the  problem  into  a  single  equation  for  the  quantity  rjf,,  given  as 
£q.  (3a)  in  the  main  text.  An  outline  of  the  derivation  of  Eqs.  (3a-c)  is  as  follows.  The 
axial  component  of  Eq.  (2c)  yields  the  constant  of  motion  ttx,«  =  7«  —  V**  —  1.  Substituting 
into  the  axial  component  of  £q.  (2a)  pves 

0^2  ~  —  V*o,,*.  (Al) 

Subtracting  this  from  Poisson’s  equation,  Eq.  (2b),  pves  Eq.  (3c)  for  p,  in  terms  of  V'*- 
Assuming  axisymmetry  {d/dO  =  0),  Eqs.  (2c)  and  (2a)  imply  us,«  =  09,,  =  0.  Substituting 
=s  +  aM,g  into  the  transverse  component  of  Eq.  (2a)  and  using  V  •  a.  =  0  and 
=  as,«  —  0  gives  Eq.  (3b)  for  ux,«  in  terms  of  Operating  on  Eq.  (2c)  with  (Vj.-) 
and  using  the  gauge  conditicm  and  the  constant  of  motion  gives 

=  Vl«„  -  ^(Vj. .  UJ.^).  (>12) 

Equation  (3a)  for  results  from  substituting  Eq.  (A2)  into  Eq.  (Al). 
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Eiquation  (3a)  can  be  written  in  the  following  form: 


The  quantities  p,  and  u,,.  can  be  written  as  functions  of  ijt,  and  |d/ using  =  y,  —rjf, — 1 
along  with  Eqs.  (2d),  (3b)  and  (3c).  Equation  (A3)  can  be  integrated  numerically  in  C 
in  a  straightforward  manner  by  starting  at  the  C  —  ^  boundary,  where  the  plasma  is 
unperturbed  and  where  =  ^$ldC  =  0  is  assumed.  For  the  purposes 

of  numerical  solution,  the  form  of  Eq.  (A3)  is  preferable  to  that  of  Eq.  (3a)  because  the 
second-order  derivative  which  occurs  in  the  last  term  on  the  right  side  of  Eq. 

(3a),  now  appears  on  the  left  side  of  Eq.  (A3).  Radial  boundary  co^^L  ons  in  Eq.  (A3)  are 
dictated  by  aadsymmetry  at  r  =:  0  and  by  the  impodtion  of  a  metallic  wall  at  r  =  >  rp. 

Equations  (1)  and  (A3),  along  with  the  auxiliary  Eqs.  (2d),  (3c)  and  (3d),  provide  a 
complete  description  of  the  laser-plasma  interaction  on  the  plasma  time  scale. 

In  numerically  integrating  Eq.  (A3)  versus  care  must  be  taken  to  avoid  a  numerical 
instability  wherein  a  suffidentiy  large  disturbance  in  if/,  with  a  suffidently  large  transverse 
wavenumber  grows  exponentially  as  a  fimction  of  This  can  be  seen  heuristically  by 
conddering  a  small  amplitude  electrostatic  wake  field  in  a  uniform  plasma  (p(°)  =  1)  with 
7a  —  1  <  <  1  and  a*, a  0.  In  this  case,  Eq.  (Al)  can  be  written 

^ - (^  +  (>M) 

where  transverse  coordinate.  For  ib^^a  <■ 

Eq.  (A4)  produces  the  usual  osdllatory  wake  field,  with  wavenumber  k^.  However,  a 
numerical  instalnlity  occurs  for  any  transverse  wavenumber  that  satisfies  >  kp/i^,{k±). 
In  jMrticular,  coodder  the  Nyquist  wavenuniber  kjy  —  w/Ax  (the  highest  wavenumber 
that  can  be  carried  on  the  numerical  grid),  where  Ax  is  the  transverse  grid  spacing.  The 
stabiBty  requirement  is  that  any  numerical  ‘Excise’'  at  k  s  ibjv  must  be  have  a  very  small 
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amplitude:  ^«(ibiv)  <  (2A*/Ap)^.  Prom  Eq.  (3c)  it  can  be  shown  the  instability  condition 
^  corresponds  to  a  negative  plasma  electron  density,  an  unphysical  situation. 

The  instability  is  therefore  numerical  in  nature  and  is  avoided  by  filtering  out  the  high 
transverse  wavenumbers  (k±  >  1/ro)  from 
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